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approach in the space-like region. We estimate their theoretical uncertainties through the use of 
the half-width rule, amounting to taking the half-width of the resonances as the deviation of their 
mass parameters. The approach embodies the meson dominance of form factors and the high- 
energy constraints from perturbative QCD. We compare the results with the available experimental 
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the experimental uncertainties, but are smaller than lattice errors. 
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I. INTRODUCTION 

Elcctroweak form factors provide valuable informa- 
tion on the internal structure of the existing composite 
hadrons, particularly on their Lorentz-invariant trans- 
verse densities [1] (for reviews of the nucleon form fac- 
tors see, e.g., Ref. [2] and references therein). On the 
fundamental level, the QCD counting rules provide the 
high-momentum behavior of the form factors [3], which 
in perturbative QCD (pQCD) acquire logarithmic correc- 
tions [4]. In recent years much attention has also been 
paid to the so-called generalized form factors, i.e., mo- 
ments of the generalized parton distributions (GPDs) (for 
reviews see e.g. [5-7] and references therein, which can 
be directly accessed on the Euclidean lattices [8], even 
when no physical experiments can. 

Electromagnetic form factors are by far best un- 
derstood, both theoretically and experimentally. An 
early and phenomenologically very successful approach 
to study these quantities was initiated by the Vector Me- 
son Dominance (VMD) model (for reviews see [9, 10]). 
On a formal level, VMD is implemented in two equiv- 
alent ways. At the field-theoretic level one postulates 
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the so-called current-field identities which state the pro- 
portionality between fields of stable mesons and the re- 
lated conserved currents with identical quantum num- 
bers. Alternately, within the framework of dispersion re- 
lations, one may saturate the matrix elements of the cur- 
rents with delta-like spectral functions. However, vector 
mesons are unstable resonances with a finite decay width, 
thus corrections to the narrow resonance approximation 
are expected. 

There arises a natural question on what numerical 
value for the mass of the resonance one should use [11- 
13]. While the rigorous definition of a resonance mass 
m R and width T R corresponds to a pole of an analyti- 
cally continued amplitude in the complex s-plane on the 
second Ricmann sheet, s = m R — iT R m R , complex ener- 
gies cannot be measured experimentally. Of course, since 
a resonance has a real mass distribution which generally 
depends on the process where the resonance is produced, 
a variety of possibilities arise to extract the maximum of 
an integrated mean value, which becomes identical and 
independent of the background in the narrow-resonance 
limit. Thus, a conservative estimate is made if the mass 
of a resonance is determined with an accuracy of about 
its half width. Following Refs. [14-18] we suggest to use 
the half-width rule: take the intrinsic uncertainty of a 
resonance mass as the range m R ± T R /2 to estimate the 
finite width corrections. The average width-to-mass ratio 
listed in the PDG was found to be (T R /m R ) = 0.12(8), 
both for mesons and baryons [16, 17]. 

From a fundamental point of view, the success of the 
VMD model remained a mystery until it was shown how 
it arises within a well defined approximation of QCD. 
Indeed, in the large-iV c limit resonances become nar- 
row [19, 20] and hadronic form factors turn out to be 
meromorphic functions under the assumption of confine- 
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ment; although in principle they contain an infinite set 
of resonances, they can generally be written as a sum of 
pole functions [21-25]. This allows us to build up an ef- 
fective theory at a purely hadronic level, with no explicit 
reference to the underlying quark-gluon dynamics: all the 
specific QCD information is contained in the chiral and 
short-distance constraints. More generally meson domi- 
nance of any vertex function with appropriate quantum 
numbers and the meromorphic property extends also to 
the case of generalized form factors. From a practical 
point of view, for the lowest-rank generalized form factors 
the only difference from the standard form factors asso- 
ciated to conserved currents is that while the momentum 
dependence remains scale independent, the normaliza- 
tion undergoes the QCD evolution [26]. Thus, general- 
ized form factors in the space-like region provide nothing 
but meson masses estimated in the large- N c limit. 

It is quite remarkable that given the simplicity of the 
lcading-iV c contribution, where only tree diagrams are 
needed, the 1/N C corrections turn out to be extremely 
complicated and despite courageous efforts [27-31] they 
have not been worked out completely. In the present pa- 
per we provide a simple way of estimating the size of the 
l/N c corrections by implementing a rather obvious idea 
that the mass of a resonance is determined with an accu- 
racy of about its half width [14-18]. While this provides 
a large- N c meaning to the half-width rule, it also yields 
rather rewarding consequences; the predicted theoretical 
uncertainties turn out to be comparable or larger than 
the corresponding experimental results, but at the same 
time smaller than current lattice estimates. In all con- 
sidered cases we find an overall agreement between the 
theory and experiment. 

The requirement of quark-hadron duality at large N c 
involves, as a matter of principle, an infinite tower of 
hadronic states [24, 25, 32-37]. This becomes clear 
for two-point functions, where further asymptotic con- 
straints between meson, the decay amplitudes and the 
meson spectra are derived . A careful scrutiny of the 
Particle Data Table confirms, with the help of the half- 
width rule [17], the radial and angular momentum Reggc 
pattern in the meson spectrum, proposed in Ref. [38]. 
This fact provides a phenomenological basis for large- N c 
Regge model calculations of form factors [14, 39-41]. 

Unlike for the two-point functions, the short-distance 
QCD constraints for the form factors may be saturated 
with a finite number of states, provided a detailed pQCD 
information (the occurrence of the logarithmic correc- 
tions) is given up. One may take advantage of this fact 
by using a sufficiently large but finite number of meson 
states, such that the correct asymptotics is reproduced 
up to (slowly varying) logarithms, which allows to im- 
pose the appropriate normalization conditions [13, 23- 



1 In fact, the only possibility to avoid the dimension-2 operators 
at large Q 2 , not present in the conventional Operator Product 
Expansion, is by assuming an infinite set of states [35, 36]. 



25, 42]. The implementation of the pQCD logs from the 
mesonic side is not at all trivial; a possible mechanism, 
involving infinitely many states, is suggested in Ref. [41], 
where it was also shown that the onset of pQCD in the 
pion form factor might possibly occur at "cosmologically" 
large momenta. We recall in this regard that almost 
model-independent upper and lower bounds on the space- 
like form factor arc established above Q 2 - 7 GeV 2 [43]. 
To be fair, it is not completely clear whether at present 
the logarithmic pQCD corrections are distinctly seen in 
the current experimental data. We note that approaches 
with an infinite number of meson resonances are consid- 
ered along the holographic framework [44, 45]. 

Usually, nucleon form factors are conveniently parame- 
terized as dipole functions, which describe the data quite 
successfully in a given range of momenta. However, this 
can only correspond exactly to a sum of two simple de- 
generate poles with opposite residues. Actually, we will 
show that if the error bars on the monopole mass are 
taken into account, one can make the dipole overlap with 
a product of monopoles within the corresponding error 
bars provided with the half-width rule. 

To summarize, our construction is based on the follow- 
ing assumptions: 

• Hadronic form factors in the space like region are 
dominated by mesonic states with the relevant 
quantum numbers. 

• The high-energy behavior is given by pQCD, and 
the number of mesons is taken to be minimal to 
satisfy these conditions. 

• Errors in the meson-dominated form factors are es- 
timated by means of the half- width rule, i.e., by 
treating resonance masses as random variables dis- 
tributed with the dispersion given by the width. 

The paper is organized as follows. In Section II we 
review for completeness the basics of meson dominance 
in the narrow width limit for the case of the nucleon. In 
Section III we digress on the role played by the finite- 
width corrections, providing a large-7V c justification for 
the intuitively obvious half-width rule for the masses. In 
Sections IV and V we carry out the analysis for several 
pion and nucleon form factors, respectively. Finally, in 
Section VI we draw our main conclusions. 



II. MESON DOMINANCE OF FORM FACTORS 

First, the implications of meson dominance on form 
factors will be illustrated with the nucleon form factors 
as an example. Quite generally, the nucleon form factors 
arc defined as the matrix element of a given current or 
composite interpolating field, J(x), 

(N(p', s')\J(Q)\N(p, s)) = u(p', s')rj(p' - p)u{p, s), 

(1) 
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where u(p, s) and u{p' , s') are Dirac spinors correspond- 
ing to the initial and final four-momentum and spin 
states, respectively. The quantity Tj(p' — p) involves 
the Dirac matrices. Lorentz indices are suppressed for 
clarity of notation. Meson dominance of the form factor 
corresponds to parameterizing the current with a super- 
position of meson fields with the same quantum numbers 
as the current, 



(2) 



which means that the meson may decay into the vacuum 
through the current, 



<0|J(0)|$ n ) =/„. 



(3) 



This also implies that the two-point correlator can be 
written as 



= E 

^— ' mi — t 



(4) 



where m n is the mass of the meson state On the 
other hand, the meson-nucleon-nucleon coupling g n is de- 
fined via 

(iV(p', s ')l(9 2 + M*)<S> n (0)\N(p,s)) = u( P ',s')g n u(p,s) 

(5) 

(g n in general involves a Dirac structure). Then 
(N(p', s')\J(0)\N(p, s)) = Yl fn(N(p', s')\$ n (0)\N(p, s)) 

n 

= u(p')F(t)u(p), (6) 
where the form factor is 

fnQn 



m = E 



mi - 1 



(7) 



It satisfies, on very general field-theoretic grounds and 
up to suitable subtractions, a dispersion relation 



1 f ImF(t') , , 



(8) 



where c.t. stands for counterterms and to is the thresh- 
old. In the narrow-resonance approximation one has the 
spectral density 



lmF(t)=TrJ2c n m 2 n 6(t-m 2 n ), 

n 

yielding, up to subtractions, the sum of monopoles 

,2 



(9) 



corresponding to Eq. (7) with c n rr? n = f n g n which is 
constant, i.e., independent oft. 



The asymptotic behavior of form factors determines 
the number of necessary subtractions. Thus, for a form 
factor falling off as <~ t~ N we have a set of conditions 

^c„m£ = k = 0,...N-l, (11) 



and thus the minimum number of meson states needed 
to satisfy these constrains is N, whence 



N 



F(t) = F(0)]J 



n=l 



mi —t 



(12) 



This simple ansatz predicts already the couplings in 
Eq. (7) "for free". One can, of course, add more reso- 
nances by multiplying Eq. (12) with a factor 



1 



dkt/ml 



l-t/r 



(13) 



where the unknown coefficient dk may be determined if 
some couplings c„ are known from the experiment. 

Quite generally, on the basis of the large-Q 2 expan- 
sion, one has [4] {-t) l+1 F t {t) - log(-i/A)"T, with the 
anomalous dimension 7^2 and weakly depending on 
the number of flavors. Fits of the form factors hardly see 
any impact of these pQCD logs at the currently available 
momenta [46]. 

The corresponding radii are given by the expansion 



F(t) 
F(0) 



1 



.V) 



(14) 



We recall that the radii are quite sensitive to chiral (1/N C - 
supprcssed) corrections and, actually, in some channels 
they diverge for m^ — > 0. 

Both two- and three-point correlators, Eqs. (4) and 
(10) respectively, require in principle an infinite number 
of mesons. Note, however, that the sign of the residues 
f n g n appearing in the form factors is arbitrary, while 
the sign appearing in the two-point correlator is positive. 
This possibly provides a quite different mechanism for 
cancellations, and hence for the form how the short dis- 
tance constraints are fulfilled. In short, the two point 
functions need infinitely many mesonic states to comply 
to pQCD, whereas the three-point functions, such as the 
form factors, can be saturated with a finite number of 
meson states. 

It is useful to notice that we may also deduce the com- 
ponent of the NN potential due to the exchange of the 
meson states <!>„, 

(N(p' 1 , S ' 1 )N(p' 2 ,s' 2 )\V\N(p 1 , Sl )N(p 2 , S2 )) - (15) 



(10) E u (Pn s 'i)9nu(pi, si)u(p 2 , s 2 )g n u(p 2 , s 2 ) 



'"7, 



t 



Via crossing, the NN scattering amplitude can be ob- 
tained as well. 
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III. FINITE WIDTH CORRECTIONS 

A question of fundamental and practical importance 
is what mass value should one use for the meson states 
in the VMD expression for the form factors [11, 17]. 
Naively, one might take the "experimental" value 2 . How- 
ever, the extended VMD formula for the form factor, 
Eq. (10), corresponds to the large- N c limit. As such, 
it is subject to the 1/N C corrections which generate a 
corresponding mass shift. The form of these corrections 
can in principle be evaluated by computing meson loops 
within the Resonance Chiral Perturbation Theory [27- 
31, 47]. The general structure of the correction for the 
form factor corresponds to the replacements 



9nfn 

ml 



t 



G n (t)fn 



rnf 



i -£(*)' 



(16) 



where E(i) is the self-energy. However, the question re- 
mains what the size of these corrections is numerically 
Strictly speaking, such a question can only be answered 
by a lattice calculation at different values of N c (see e.g. 
Ref. [48] and references therein). Unfortunately, as we 
argue below, within a purely hadronic resonance theory 
we can only make an educated guess, since there are un- 
determined counterterms encoding the effects of the high- 
energy states not considered explicitly. For instance, for 
the case of the p-meson we may take into account the de- 
cay into 2tt, which is a real process, but also the virtual 
KK excitation, etc. Our lack of an explicit knowledge 
on all excitations makes it difficult, if not impossible, to 
predict the mass shift reliably. 



A. Mass shift and the width 

To elaborate on the mass-shift effect in a greater de- 
tail, let us consider the two-point function yielding the 
mesonic resonance propagator, 



D(s) = 



1 



E(«) 



(17) 



The mass parameter m is the tree level resonance mass, 
which is 0(N®), whereas the self energy, coming from 
meson loops, is suppressed, £(s) = 0(N < 7 1 ). 

Let us consider, for instance, the self-energy correction 
of the scalar or vector mesons due to pion loops. Ana- 
lyticity implies that the self-energy satisfies a dispersion 
relation 3 



S(s) = c.t. + - f 

7T J 4 



* d /^' +i ° + \ (18) 

Ami S S 



where c.t. means suitable subtractions. The pole posi- 



tion, s — sr ; = m R — imuTji, is given by 
sr - ml - S(sij) = 0. 



(19) 



This is a complicated self-consistent equation, but within 
the 1 /N c expansion it can be solved perturbatively, yield- 
ing 

SR = m 2 + 2m Q AmR-iTRm + O(N- 2 ), (20) 



where 



Amp = 



2too 



ReS(m^), 



T R = Im£(m§). 

m 



(21) 
(22) 



Note that the imaginary part is proportional to the cor- 
responding decay width, 

In general, there appears a threshold momentum de- 
pendence for the decay amplitude which is proportional 
to the phase space. The form reflects the spin of the 
resonance, such as 



T(s) 



P(s) 



p{m R ) 



2J+1 



(23) 



with p(s) = \fl — 4m 2 / s = p/\fs, where p is the center- 
of-mass momentum when m§ — > sq. Obviously, the num- 
ber of subtractions in Eq. (18) depends on the assumed 
off-shcllness. For the previous choice of T(s), which be- 
comes a constant at large s, we need at least one sub- 
traction, which we may choose to be, e.g., at s = and 
thus in terms of the principal value integral we have 



Ami 



2m 



i r 

ReS(O) + — I ds 

n Jim! 



,m\ ImS(s') 



(24) 



Therefore AmR depends on an arbitrary constant 
ReS(O) = 0{N^ 1 ), which cannot be determined from the 
dispersion integral or the lowest-order parameters and 
hence naively becomes independent of the width. Other 
momentum-dependent widths, not vanishing at high s, 
may introduce additional subtractions. The present dis- 
cussion illustrates our statement that one cannot gener- 
ically compute the mass shift in a model-independent 
way 4 . 

This lack of predictive power within the purely 
hadronic theory is not surprising. However, from a mi- 
croscopic point of view the meson self-energy can be un- 
derstood as the coupling of the qq bound state to the 
meson continuum and physical resonances turn into Fes- 
chbach resonances. The relevant scale corresponds to the 
string breaking distance, defining a physical momentum 
scale which may be described as a transition form factor 
from qq states to mesonic channels. This implies that the 



2 This value also depends on the experimental process and may 
differ within the half-width rule. 

3 We disregard spin complications, see, e.g., [49] for details. 



4 This is so provided no further information is available. 
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mass shift due to closed channels is necessarily negative 
as it corresponds to second-order perturbation theory be- 
low the closed channels, but also that the mass shift due 
to the open channel, scales exactly as the decay width. 
In Appendix A we analyze some specific models where 
we can see that within uncertainties a natural rough es- 
timate of the mass-shift is given by the half-width rule. 



B. Finite width effects in the space-like region 

Finite width corrections for the pion charge form factor 
were pioneered by Gounaris and Sakurai [51]. They have 
implemented the e + e~ — > tt + tt~ hnal-state interactions 
in the timelike region, where they are crucial. In this 
section we analyze the influence of widths on the space- 
like region. To this end, we use Watson's theorem on 
hnal states, which can be written as 5 



F y (s + iO+) _ Tii(s-iO+) 
F y (s-iO+) ~ T n (« + iO+) 



e 2i«ii(«) ) 4m 2 < s < 4m |_ ; 



(25) 



where Fy(s) is the form factor, Tn(s ± i0 + ) is the tttt 
partial-wave scattering amplitude in the vector-isovector 
channel ( J, I) = (1,1), and 5n(s) is the corresponding 
phase-shift. A well known solution to this discontinuity 
equation is given in terms of the Omnes function, 



0(,s) = cxp 



i r 

n J 4m 



Sj(s') 



2 s'(s-s') 



(26) 



which fulfills O(0) = 1. A solution to Eq. (25) is given 
by just taking Fv(s) = P(s)Cl(s), with P(s) being an 
arbitrary polynomial. Choosing P(s) = 1 we have 



F v (s) = Q(s). 



(27) 



In the case of a zero-width resonance the pahse-shift is 
£li(s) = tt/2 6(s — ra 2 ) and one gets the monopole form 
factor, 



Fy( S ) 



mi 



mt> — s 



(28) 



featuring VMD in its simplest version. 



5 Another way of writing the relation, which generalizes trivially 
to coupled channels, is through the use of the Bethe-Salpeter 
equation FT" 1 = TV [52], which yields DiscF v (s + if^)" 1 = 
DiscTn(s + iO+)- 1 . Then 

DiscF; (s + i0+ ) " 1 = Disc ^2 T »i ( s + i0+ ) ~ 1 ■ 

3 

Above the KK production threshold Watson's theorem requires 
considering also the Kaon form factor and the corresponding 
extension to coupled channels, i.e., the mixing with tttt — > KK 
transitions. 



We use a simple Breit-Wigner (BW) parameterization 
for the vector-isovector 7T7r-phase shift, obtained as 



2M llW = Dy(s + lO + ) 

D v (s-iO+) ' 



where 



[D v (s)} 1 = s-m 2 p + im p T p 



(s - AmDrrip 



Am 2 )t 



(29) 



.(30) 



The p-wave character of the p — > 2ir decay can be rec- 
ognized in the phase space factor. The Omnes form fac- 
tor is depicted in Fig. 1. As we can see, the finite width 
correction lies within the band corresponding to the half- 
width rule imposed on top of the monopole form factor, 
considering here m p = 0.77 GeV and T p = 0.15 GeV. 



C. The half-width rule 

As we see, the sublcading 1/N C corrections in the 
space-like region essentially correspond to keeping the 
meson dominance form and to changing the parameters. 
By making simple calculations we have seen that a con- 
servative bound on the mass shift is given by the half- 
width rule. From a spectral point of view this is quite 
natural when we appeal to the Kallen-Lehmann repre- 
sentation of the resonance two-point function, 

™-r*l^w (31) 

Wc may then use the probabilistic interpretation of the 
line shape 



P(p) = Zp{n). 



(32) 



If we take a Breit-Wigner shape for D(s) (neglecting the 
threshold effects), we have 



-Pbw(m) 



2r M 2 



7T Ou 2 - m 2 ) 2 + rv 



(33) 



which is normalized to unity, / dpP(p) = 1. 

The random implementation for a given distribution 
is obtained in a standard way by inverting the relation 
expressing the coordinate independence of probabilities, 



P(/j,)dp = dz, 



(34) 



with z € U[0, 1] being a uniformly distributed variable 6 . 
The result for a BW shape in the case of the p-meson for 
N = 10 4 samples is shown in Fig. 2. The idea amounts 
to treating the resonance mass as a random variable and 



6 For a Lorcntz distribution, P(/x) = Z^/{(p? - M|) 2 + r 2 R M%), 
the relation is given by fi 2 = + M^Vr tan (ttz/2), with 



2 6 U[0, 1]. This distribution must be cut and normalized when 
negative ft 2 values are generated. 
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0.0 0.5 1.0 1.5 2.0 ' o 1 2 3 4 

V7 [GeV] Q 2 [GeV 2 ] 

FIG. 1. (color online). The (7, J) = (1,1) 7T7T scattering phase shift as a function of the center-of-mass energy (in GeV). We 
use the BW representation discussed in the text. The data are from the analysis of Ref. [50] (left panel). The pion charge 
form factor in the space-like region as a function of the momentum Q 2 (in GeV 2 ) for the Omnes representation (solid red), 
compared with the simple monopole form (dashed blue). The band corresponds to taking a monopole with the half- width rule 
(right panel). 




0.0 0.5 1.0 1.5 2.0 ' o.O 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

m p [GeV] q 2 [GeV 2 ] 

FIG. 2. (color online). Sampling of the p-meson mass according to the BW spectral distribution. We sample N = 10 4 values 
and bin them with Am = 20 MeV (left panel). The monopole form factor F 7I7I {Q 2 ) = m 2 /(rn 2 + Q 2 ) sampled with the previous 
distribution and compared to the experimental data [53-63] (right panel). 



to propagating its effect in all obscrvables. Of course, 
different shapes produce somewhat different confidence 
levels. For definiteness, we will take Gaussians which 
have shorter tails and are symmetric around the reso- 
nance value. 

In Fig. 2 we plot the monopole form factor 



F V (Q 2 



(35) 



mj + Q 2 

according to the BW distribution of the mass. 

IV. PION FORM FACTORS 
A. Electromagnetic form factor 

The charge form factor of the pion is given by 

{n + (j/)\JJT(0)W + (p)) = (p^+P^Mq 2 ), (36) 

with q = p' —p and the electromagnetic current J° m (x) — 
J2 q=u ,d,s,... e qQ( x )lfj.l{x), where e q denotes the quark 



charges in units of the elementary charge. The charge 
normalization requires 



F V {0) = 1. 



(37) 



Actually, in the space-like region, where t = —Q 2 , F(t) 
is real and at large Q 2 values the pQCD methods can be 
applied, yielding asymptotically [68-73] 



Fv{-Q 2 



16irf 2 a(Q 2 ) 



Q 2 > m 2 



1 + 6.58 



a(Q 2 



(38) 



with f n — 92.3 MeV denoting the pion weak decay 
constant, and m standing for the lowest vector meson 
mass. If we ignore the slowly varying logarithm, we get 
Fy(t) = 0(i _1 ) and in the large N c limit one has 



Fv(t)= J2 

V=p,p', 



c v — T 



(39) 



where J^v Cv = 1 an( ^ c v ~ 9VnnFv /my ■ 
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FIG. 3. (color online). Top row: The pion charge form factor with the half-width-rule compared with the experimental data 
[53-63] (left) and the pion-photon transition form factor compared with the experimental data [64-67] (right). Bottom row: 
The pion charge form factor multiplied by Q 2 (left) and the pion-photon transition form factor multiplied by Q 2 (right); the 
horizontal dashed line represents the asymptotic value 2f 7r - 
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FIG. 4. (color online). Axial form factor in the space-like region t < 



The simplest formula fulfilling this constraint and (37) 
is the VMD solution 



F v (t) 



m 2 — t ' 



(40) 



whence g p7rK f p = m p . The p — 7 coupling is given by 



whereas the p — > mr decay is 

r(p — > 7T7r) = 



3 m p ' 



9 P ^rn p 



48tt 



(41) 



(42) 



Using the PDG numbers one gets for m p = 0.77 GeV the 
values f p = 0.156 GeV and g p7V1T — 6. 

One of the important features of the pion form factor is 
that the radius has large chiral corrections, thus we may 
improve the phenomenology by adding one extra meson, 
p'. Then 



F v (t) = (l-c)-^- 
mi 



ra 2 , 



2 

m 2 , 



such that 



h r 2 ) = (l-c)^+c ] 



m 2 , 



(43) 



(44) 



Thus imposing the physical value of (r 2 ) we get the 
value of c for any m p and m p > . Taking again the PDG 
values for those quantities (m p = 0.77549(34) GeV, 



= 1.465(25) GeV, 



(0.672(8)fm) 2 ) we obtain 



c= -0.227(39). 

One could also carry out the analysis the other way 
around, starting from Eqs. (40,43) with all the constants 
treated as free parameters to be determined by a fit 
to the experimental data. That way it was shown in 
Refs. [11, 37] that one can retain a precise value for (r 2 ), 
even though the masses do not have precisely their phys- 
ical values. In the large- N c limit, when the functions be- 
come meromorphic, this fitting procedure is mathemat- 
ically safe thanks to the convergence theorems from the 
Pade Theory [37]. In this framework, Fig. 3 of Ref. [11] 
shows the half-width-rule as a good estimation of the 
systematic error done on the determination of the poles 
when fitting the space-like data. 



B. Axial form factor 

The axial form factor of the pion is intimately related 
to the pion radiative decay — » I 71/ (with I stand- 
ing for e or fj,) and its hadronic contribution. The decay 
proceeds via ordinary inner bremsstrahlung (IB) from 
the weak decay — > l^v accompanied by the pho- 
ton radiated from the external charged particles, and the 
structure-dependent interaction (SD) between the pho- 
ton and the virtual hadronic states, with contributions 
of both vector and axial-vector form factors. 

For the transition ir + —¥ "fv e e + [74], the structure de- 
pendent amplitude is given by 

M SD = ieG F cos 8 < c u ulp {l - T5 )u e e*M^/V2m 7r ,(45) 

where e M is the polarization vector of the photon, Gf 
is the weak interaction coupling constant, and 8 C is the 
Cabibbo angle. The hadronic contribution is enclosed in 
the amplitude M^ v and reads: 

M"" = F v {t)e^p a q p 

-iF A (t)[q v (q^+p^) - <Tq ■ (q + p)\ , (46) 

with Fy :A (t) denoting the vector and axial- vector form 
factors, respectively. By assuming axial-meson domi- 
nance, we have 



M 2 

F A (t) = F A (0) A 



(47) 



with F A (0) = 0.0119(1) [75]. The result obtained with 
the half-width-rule is presented in Fig. 4. 



C. Transition form factor 

The pion-photon transition form factor 7r° — > 77* has 
been subjected to vigorous discussion in recent years. 
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FIG. 5. (color online). Band: the pion-photon transition form 
factor of Eq. (52). Points: various experimental data [64-67]. 
The horizontal line represents the theoretic asymptotic value 
of 2/,. 



Firstly, its value at the origin is fixed by the chiral 
anomaly, 



F(0) 



1 



47T 2 /,' 

while its asymptotic behavior is given by 



F(Q 2 ) -> 



6/ ff 

NrQ 



A simple model fulfilling both conditions is 

2 1 m p 

F{Q ) = 4 7 r 2 / 7r m 2 + Q 2 ' 



provided one has the relation 



247T 2 / 2 
Nr ' 



(48) 



(49) 



(50) 



(51) 



which gives m p = 823 MeV for /„■ = 92.6 MeV or m p 
770 MeV for f n = 86.6 MeV in the chiral limit. 

If we include two resonances [76], p and p' , we get, 
after imposing the anomaly and large-Q 2 behavior, 



F(Q 2 



m 2 p m 2 p , 



-24/> 2 Q 2 /7V ( 



^ 2 U 



(m 2 



Q 2 )(m 2 ,+Q 2 ) 



(52) 



The result is shown in Fig. 5, using m p = 0.775 GeV, 
m' p = 1.465 GeV, T p = 0.150 GeV and T' p = 0.400 GeV. 

One could even go beyond this approximation includ- 
ing a third resonance (the p"). That introduces a new 
parameter that can be fixed by the derivative of the form 
factor at the origin, the parameter [13]: 



F(Q 2 



(53) 



99 9 
1 m p m p ,m 



p" 



bQ 2 + 24/ 2 7r 2 Q 4 /^c 



47T 2 / T (TO 2 



2 )K 



Q 2 )(m 2 „ 



+ Q 2 ) 



where the parameter b can be obtained through a match- 
ing procedure to the low-energy expansion of F(Q 2 ), 
i.e., b = TO 2 m 2 ,TO 2 „(-% + A- + At H a-)- Given 
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= 0.135 GeV, m p ,m' p , m!' p = 1.720(20) GeV, and 
= 0.032(4), we obtain b = 5.82(18). 
Figure 3 of Ref. [13] shows how the half-width-rule 
provides a good estimate of the systematic error on the 
determination of poles of rational approximants, such as 
Eqs. (52,54), when fitting to the space-like data [64-67]. 



D. Gravitational form factor 

The gravitational quark form factors of the pion [77] , 
0i and 2 , are defined through the matrix element of 
the quark part of the energy-momentum tensor in the 
one-pion state, 



(7T b ( P ') | 0^(0) | 7T») = 

V [(g^q 2 - qXqne^q 2 ) + 4P^P»Q 2 (q 2 



(54) 



where P = \{p' + p), q = p' — p 1 and a, b are the isospin 
indices. The gravitational form factors satisfy the low- 
energy theorem ©i(0) - 2 (O) = 0(m\) [77]. The trace 
part is dominated by scalar states, while the traceless 
component is dominated by spin-2 tensor mesons. 

Following the standard notation for the moments of 
the pion generalized parton distributions (GPDs) [26], 
we introduce: 



1 



A 20 {t) = Mt) 



A 22 {t) 



1 



2 (t), (55) 



where the symbols 0,(i) denote the quark parts of the 
gravitational form factors of Eq. (54). We consider also 
the moment 



A 10 {t) = F v {t) 



(56) 



with F v (t) denoting the electromagnetic form factor de- 
scribed in Sec. IV A. 

The results of monopole representations with the half- 
width rule are shown in Fig. 6, where we consider 



A ^{t) = ^ptl^ ^o(t)=A 20 (0)^- t , (57) 

P ' h 

with m h = 1.320 GeV, T h = 0.185 GeV, and A 20 (0) = 
0.261. We also show in Fig. 6 the low-energy theorem 
that relates A 2 o(t) with the other moment, A 22 (t). The 
relation A 22 {§) = —jA 20 (0) is compared to the lattice 
data of Refs. [78, 79]. 



V. NUCLEON FORM FACTORS 
A. Electromagnetic form factors 

In the non-strange sector the electromagnetic current 
is given by 



1 



(58) 



where Jg(x) is the baryon current and J v 3 (x) is the third 
component of the isospin current 



J£(a0 = ?»f?(i), 

The matrix elements of these currents are 
(N(p')\r B (0)\N(p)) = 



(59) 



u{p') 



rF(=°( q 



2M N 2 U ' 



*(p)> 



(N(p')\J v a (0)\N(p)} 



u{ P 'y 



w^q« p i=i 
2M, 



N 



if) 



(60) 



where q = p' — p, and F± and F 2 are the Dirac and Pauli 
form factors, respectively. The relation to the proton and 
neutron form factors is 



Ff 
F? 



( F I=0 

(Ft 



Ft 1 ), 
Ft 1 ), 



where 



Ff (0) = 1, 
ff(0) = k p , 



F?(0) = 0, 
F 2 "(0) = K n . 



(61) 



(62) 
(63) 



The quantities n p = 1.793 and k„ = —1.913 are the 
anomalous proton and neutron magnetic moments, re- 
spectively. The electric and magnetic Sachs form factors 
are defined as 



G P E (q 2 ) = F?(q 2 ) + 



Gl(q 2 ) 



4M 2 N 



Fi(q 2 ), 



Fl{q 2 )+F^q 2 l 
G n M (l 2 ) = F?(q 2 ) + FZ(q 2 ). 



G P M (q 2 ) 



The normalization conditions become 

G*(0) = 1, G n E (0) = 0, 
G p M (0) = p p , G n M (0) = »n, 



(64) 



(65) 
(66) 



where p p — 2.79pN and jj, n = — 1.91/Ujv, with p^ — 
e/(2Mjv) denoting the nuclear magneton. 

The asymptotic behavior for t — > — oo is given by [3] 

f+'Fiit)^ [logH/A 2 )] -7 , (i = l,2) (67) 

where the anomalous dimension 7 <~ 2 is weakly depend- 
ing on the number of flavors. As mentioned before, such a 
slowly changing log behavior cannot be reproduced with 
a finite number of resonances, thus we assume it to be 
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FIG. 6. (color online). Electromagnetic form factor Aio(t), the corresponding (—t)Aio(t), and the quark part of the gravitational 
form factor A2o(t) and (— i)A2o(t), compared to the lattice data from [78] (top and middle row respectively). Spin-0 gravitational 
form factor of the pion, A22 (t) , from the lattice calculation of Refs. [78, 79] extrapolated to the physical pion mass. Red squares 
correspond to A22(t) and blue circles to — Aw(t)/A (bottom row). 



constant. At the leading order in the large- N c expansion 
the form factors read 



gujNNfui-y 



^ J TTi — 7 



V 

V P 



g P NNf P 

m 2 p — t 

fpNNffry 



v m l- t 



(68) 



We define the strong isoscalar and isovector vertices 



(N{p')\^\N{p))=u{p') 



(N( P 'M\N(p))=u(p')^ 



9uiNNl^ + fuiNN 
gpNNl^ + fpNN 



2M N 



2M 



N 



u(p). 
(69) 



According to Eq. (67), the minimum number of reso- 
nances is two and three for the Dirac and Pauli form 
factors, respectively. We use the normalization condi- 
tions at the origin, the asymptotic conditions, and fix 
the vector and tensor couplings to the lowest-lying res- 
onance, Qujnn 1 fuNN, 9 p nn, fpNN m the isoscalar and 
isovector channels. Our illustrative goal here is to predict 
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the form factors without attempting a detailed fit to the 
data. 

After all the conditions are imposed, we get 



F{=\t) = 

Ft°(t) = 
Ft\t) = 



1 



1 - c t/ml 



2{l-t/ml){l-t/ml,){l-t/ml„y 



1 - cxt/m 2 ,, 



2(l-t/m 2 p )(l-t/m 2 p ,)(l-t/m 2 p „y 
1 1 
2(l-t/ml)(l-t/ml,)(l-t/ml„y 
1 1 
2(l-t/m 2 p )(l-t/ m 2 p ,)(l-t/m 2 pll )> 



(70) 



where the constants cq and c\ are determined from the 
values of g U NN an d 9 p nn, respectively, as 

g^NNL-f _ 1 1 - c m 2 /m 2 ,, 

2(l-m 2 /m 2 ,)(l-m 2 /m 2 ,,)' 
1 l-c m 2 /m 2 „ 



9pNNfp 



2 (1 - m 2 p /m 2 p ,)(l - m 2 /m 2 „) ' 



(71) 



Thus, with the electromagnetic decay widths for the vec- 
tor mesons, F p = 0.149 GeV, we find (Eq.(41)) T(p -> 
e + e _ ) = Aira 2 F 2 /m p /3, which is numerically equal to 
6.4 KeV. 

Detailed fits implementing VMD [80] require large OZI 
violations as well as huge departure from the flavor SU (3) 
symmetry (see also [81]). In particular, the value of 

9ujNN ^ 9^nn [46, 82]. However, the simple Kelly pa- 
rameterization [83] provides a successful fit in the space- 
like region as a rational function with the correct large 
momentum behavior. 

The result of varying the masses according to the half- 
width rule is presented in Fig. 7. We consider the SU(3), 
case where g U NN / 9pNN = fuNN / fpNN = 3, as well 
as a 30% violation for the ratio (orange and blue band 
in Fig. 7, respectively). We recall that in the meson- 
exchange models [84] or in dispersive analyses of the nu- 
clcon form factors [46, 82] even much larger symmetry 
breaking is needed to comply with the phenomenology. 



B. Axial and pseudoscalar form factors 

The axial matrix element of nucleon is defined by (for 
a recent review see, e.g., Ref. [88] and references therein) 



(N(p')\J^(0)\N( P )) = 



«(p')y75 



u(p), (72) 



where Ga(o 2 ) and Gp(q 2 ) are the axial and the induced 
pseudoscalar form factors, respectively. The QCD axial 
current is 



(73) 



while current-field identity relating it to the axial-meson 
field A^ v and the pseudoscalar field P is 

JaW = E fAd v A a ^{x) + fpd^P a (x). (74) 

Therefore we have 

A A 

„ sr^ ^ m nIa9ann , AM N F P g PNN 



, M\-t 



M 2 -t 



(75) 



We use here the extended PCAC form [89] , which for the 
on-shell mesons reads 

«>) = E.^ M p p >)' ( 76 ) 

yielding the following relation among the form factors: 

2M N G A (t) + ^-GpW = £ 2 -^P^9pnn. (77) 

The pseudoscalar-nucleon coupling is defined by 

(p'\(d 2 + M 2 )P»{x)\p) = g PNN u{p>)i lb T a u{p). (78) 

/From here we get the (extended) Goldberger-Treiman 
relation 

MpjgA — £ F P g PNN = f^g-nNN + fn'9n'NN + 

P 

(79) 

The high-energy behavior of the weak form factors in 
QCD was discussed many years ago [90, 91]. At high 
Q 2 , one has for the isovector [92] and isoscalar [93] the 
asymptotic behavior, 



Q A G A {-Q 2 ) -> const. 
We also have the sum rules 

9A = £ fAgANN, 
A 

= YjfA9ANNM\. 



(80) 



(81) 



It is noteworthy that most determinations of the axial 
form factor proceed via a dipolc fit, 



G A (t) 



9A 



(l-t/Ai) 2 ' 



(82) 



suggesting a 1/t 2 fall off at large t. The values of 
the parameter are Aa = 1.026(21) GeV or A^ = 
1.069(16) GeV, depending on the process [88]. In the lit- 
erature Aa is denoted and called the axial mass (see, e.g., 
[88]). This is not our axial meson mass, since Eq. (82), 
although phenomenologically successful, cannot be justi- 
fied from a field-theoretic point of view and is in contra- 
diction with the large- N c motivated parameterization. 
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FIG. 7. (color online). The electromagnetic nucleon form factors compared to the data for the proton [85] and neutron 
(and references therein). Orange band: q^nn = 9, blue band: Quinn = 12. 
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FIG. 8. (color online). The isovector-axial nucleon form factors using two axial masses (left panel). The isovector-pseudoscalar 
nucleon form factor using two axial masses and the lightest pion considered in [87] (right panel). 



The minimum meson-dominance ansatz compatible 
with low- and high-energy constraints reads 



G A (t) = g a 



9 ? 

m„ m i 



G P {t) = G A {t) 



Gp(0) 



ty 



— t 



(83) 



By applying the half-width rule to this parameteriza- 
tion, i.e., using m ai = 1.230 GeV, m a i — 1.647, r ai = 
0.425 GeV, and = 0.254 GeV, we get the results 
depicted in Fig. 8. As we can see, the results are in rea- 
sonable agreement with the data. Actually, the two axial 
mesons are incorporated as a product of monopoles, but 



since they have an overlapping spectrum, the net effect 
is essentially a dipole form factor with an average mass 
which is somewhat larger than the usual dipole cut-off. 



C. Gravitational form factors 

The discussion of the nucleon gravitational form fac- 
tors follows closely the pion case with suitable changes. 
For the nucleon case, the quark contributions to these 
form factors have been determined by the QCDSF Col- 
laboration [94] and the LHPC Collaboration [95]. 

The decomposition, corresponding to the energy- 
momentum tensor matrix elements taken between nu- 
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cleon states, reads 

(p'|e«» = u(p' 

+ BUt) 



CUt) 



A„A 



v2l 



A' 



«(p), (84) 



where er^" = |[7 M ,7 1/ ] (the Bjorken-Drell notation), the 
scalar functions are moments of the GPDs, the momen- 
tum transfer is denoted as A = p' — p, and the average 
nucleon momentum is P — (p' +p)/2. Taking the trace 
and applying the Gordon identity, 2MNu(p')"f IJ 'u(p) = 
u(p')(i<r ,ip A p + 2PP)u(p), as well as the Dirac equation, 
(j$ - M N )u{p) = and u(p')(]f - M N ) = 0, we obtain 
the following expression for the spin-0 gravitational form 
factor of the nucleon: 



e%(t) = m 



N 



3f 



A I 



T^2o(^) 
A 



(85) 



whereas the spin-2 (normalized) component becomes 

AUt) 



F«{t) 



AIM 



(86) 
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FIG. 9. (color online). Spin-0 gravitational form factor of 
the nucleon, Ge{t) obtained from the lattice simulations of 
Ref. [95] at the pion masses m w — 352.3MeV (blue circles) 
and m-n = 356.6 MeV (red squares). 



The trace is dominated by scalar, ++ , states. For the 
monopole representation of the gravitational form factor, 



G e (t) 



Jo 



fo ~ t 



(87) 



depicted in Fig. 9, we use the updated values of the sigma 
meson properties from the latest edition of the PDG Ta- 
bles [75], i.e., m a = 475 MeV and T a = 550 MeV. 

The traceless part of the EM tensor corresponds to 
a spin-2 isoscalar gravitational form factor which natu- 
rally couples to the ji meson within a tensor dominance 
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FIG. 10. (color online). Spin-2 isoscalar gravitational form 
factor of the nucleon, Fr(t) obtained from the lattice simu- 
lations of Ref. [95] at the pion masses m n — 352.3 MeV and 
m-n- = 356.6 MeV (blue circles and red squares respectively), 
together with the results of Ref. [95] linearly extrapolated to 
the physical pion mass (black triangles). 



approximation. For the nucleon this FF has been deter- 
mined by two lattice groups in Rcfs. [94, 95]. Actually, 
in Ref. [94] a dipole fit describes the data successfully, 
namely 



F T {t) 



1 



with At = 1.1(2) GeV, if a linear extrapolation in to 
the physical point is assumed. Assuming an asymptotic 
falloff for the form factor, such that Fy(t) = 0{t~ 2 ), we 
just take a sum of two monopoles that reduces to 



F T (t) 



>■ r 

J-2 



mi — t m 2 t , — t 
J2 f? 2 



(89) 



The PDG Tables quote m h = 1.320 GeV and T h = 
0.185 GeV, and for the first excited state gives mp = 
1.525 GeV and T fi = 0.073 GeV. 7 

As we can see from Fig. 10, similarly to the case of 
the axial nucleon FF, the dipole FF with an uncertainty 
essentially corresponds to two monopoles after the half- 
width rule has been implemented. 



VI. CONCLUSIONS 

In the present work we have taken advantage of the 
well-known fact that in the large-7V c limit of QCD the 
generalized hadronic form factors, probing bilinear qq op- 
erators with given J PC quantum numbers, feature gener- 



7 This is somewhat different from the results of Ref. [96] with about 
four times larger width for the /2(1565) meson. In any case, for 
the spin-2 form factor the difference becomes rather irrelevant in 
the region below 1 GeV. 
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alized meson dominance of qq states with the same quan- 
tum numbers. They assume the monopole form, 

AB 2 

(A(p')\J(0)\(B(p)) ~Y1 S^T' ( 9 °) 
z — ' mt — t 

n n 

where m n are the meson masses and c^ B suitable cou- 
plings. Thus generalized form factors at some finite mo- 
mentum transfer essentially measure the masses of the 
lowest lying mesons. 

The goal of this paper was to present a comprehensive 
analysis of the pion and nucleon form factors, providing 
a theoretical uncertainty following from the half-width 
rule. We have incorporated 

• The correct asymptotic power like behavior of form 
factors (short-distance constraints). 

• Low-momentum normalization constraints. 

• Minimum number of mesons with the relevant 
quantum numbers in each channel. 

• Theoretical error estimates based the half-width 
rule. 

Given the approximate nature of the underlying large- N c 
expansion, we should not expect perfect agreement with 
data. Rather, the addressed question is how we can esti- 
mate the accuracy of the large- N c expansion for hadronic 
form factors, which are basic experimental quantities. In 
the present paper we provide arguments in favor of the 
simple rule, where a rough estimate is given by vary- 
ing the resonance mass within its width range. As we 
have seen, this provides a surprisingly close answer to 
the data, which fall within the bands produced with the 
half-width rule. While this presumably is a conservative 
assumption, it still provides a cheap estimate based on 
an independent source of information. 

Extension to other baryons and mesons is straightfor- 
ward, as well as to the transition matrix elements. One 
useful application of the meson-dominance scheme is the 
a priori determination of Generalized Form Factors, just 
based on the PDG Tables, for which abundant data start 
to be produced on the lattice. Our calculations show that 
improving on the uncertainty predictions based of the 
half-width rule may require highly refined lattice studies 
beyond the present accuracy. 
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FIG. 11. (color online). Mass-shifts Am (in GeV) of the p- 
meson (top panel) and the a- meson (bottom panel), due to 
pion loops, plotted as functions of the cut-off A (in GeV) for 
several cut-off functions: monopole (dashed brown), dipole 
(solid red), sharp cut-off (dot-dashed blue). 



where pn = p(m R ) = y/l — im 2 /m R . Note the addi- 
tional s factors appearing in the widths. By using the 
dispersion relation for the self-energy, Eq. (18), we can 
see that the real part of the integral is quartically di- 
vergent and thus three subtractions are needed. This is 
equivalent to fixing the mass and the width. On the con- 
trary, the imaginary part is finite and provides the width 
at the physical value. This argument shows that we need 
more input information and until then have no predic- 
tive power. Within a large N c environment [97] we use 
BW and not pole reference subtraction values for both 
mesons in this study. 

In order to get an estimate, we assume some form of 
the hadronic form factor for the vertex pnir. Out of ig- 
norance, we consider different cut-off functions 



Appendix A: Cut-off dependence of the mass shift 

In this appendix we analyze the mass shift of the low- 
est lying scalar and vector resonances due to the most 
relevant decays S — > tttt and V — > tttt. With a chiral 
derivative coupling, one obtains [97] 



r*(s) 



3ms s 

167T/4 



PS 



Cd + (Cm - c d) 



2mi 



G^mys 3 



(Al) 



GpT T7: (q 2 ) 



A 5 



s + A 2 



(A2) 



which preserve the imaginary part, G pn7r (m 2 ), and hence 
the width. Here n — 1 corresponds to a monopole, n = 2 
corresponds to a dipole, and the limit n — > oo corre- 
sponds to a sharp cut-off. We naturally expect A in the 
range above m p and around m p >, which corresponds to 
ignoring all states above these energies. With these con- 
ditions, and assuming a small correction, we get a mass 
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shift 

Am p = 1 I" ds G P^Ej>^A . (A3) 
7T J 4m 2 s{s - m 2 p ) 

The results are depicted in Fig. 11 as a function of the 
cut-off scale. As expected, the overall order of magnitude 
is quite compatible with the half-width rule. 



In a theory with a finite number of resonances there is, 
of course, an implicit high energy cut-off corresponding 
to the next (not included) meson. Assuming a vanishing 
contribution of the higher energy states, one may make 
a rough estimate of the mass shift. For instance, for the 
p-meson one has [98] a shift of about half the width. Of 
course, this may be partly a numerical coincidence, but 
it illustrates the point that parametrically the mass-shift 
and the width scale in a similar way. 
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